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Abstract: DeshDande & Rankin (1999, 2001) claim that the frequency of the very 



narrow feature, in the spectrum of radio flux variations of PSR B0943+10, is an 
alias of its actual value. For this, they need to interpret the narrow feature as a 
very stable pattern of drifting sub pulses on t he p olar cap of this pulsar, within the 



framework of the Ruderman & Sutherland (1975) model. This paper discusses 



several factors within the same framework, that may perturb a very steady pattern 
of drifting sub pulses, leading to a broadening of the spectral feature. One should 
be cautious in (1) identifying such very narrow spectral features with a very sta- 
ble pattern of drifting sub pulses, and (2) use this picture to resolve the issue of 
aliasing in PSR B0943 + 10. 

Keywords: pulsars: general - pulsars: individual (PSR B0943+ 10): stars - neu- 
tron — fluctuation spectrum — drifting sub pulses. 

1 Introduction 

The rotation powered radio pulsar PSR B0943 + 10 has a very narrow feature in 
the longitude resolved spectrum of intensity fluctuatio ns; its Q, defined as its cen 
tral frequency divided by its width, is relatively high (TTavlor & Huegueninlll971 



tral frequency divided by its widtn, is relatively nign ( laylor & hlueguenin ly /l; 
Backer. Rankin & Campbell [l97.4lsieher & Osteill 975h. RecenflvlPeshpande & Rankin 
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( 19991 1200 lb henceforth DR1999 and DR2001) put its Q at » 500, and also 
claim that this very narrow spectral feature, occurring at 0. 465 cyc l es p er pulsar 
period [cpp], is an alias of the actual value 0.535 [cpp]. iBackei dl973h noted 
such a general possibility in radio pulsars and claimed, based on the phase infor- 
mation in the fluctuation spectr um, that PSR B2303+ 30 has an aliased spectral 
feature. This was contested by Sie ber & Ostej dl975l) . who st ate on their p age 
326 that "even a phase analysis, contrary to Backer's statements (Backe ll973h . is 
unable to decide between the two possibilities", viz., whether the spectral feature 
is aliased or not. Indeed, in a later paper Backer did not repeat such a claim for 
PSR B0943+10, whose fluctuation spectrum is similar to that of PSR B2303+30 
(backer. Rankin & Campbelll! 1975b: he quoted the the t rue value as 0.465 [cpp]. 
However, In the view of DR20Ql. lSieber & Ostei ( 1975 ) "came to the wrong con- 
clusion". 

Like earlier workers, DR2001 associate this high Q spectral feature with the 
well known phenomenon of systematically drifting sub pulses in radio pulsars, and 
offer three arguments, involving signal processing of position modulated pulses, 
which are by themselves not sufficient t o resolve the aliasing issue. Th ey need 
to make the crucial assumption that the iRuderman & Sutherland ( 1975b model 
(henceforth RS1975) applies to the data of PSR B0943+10. Only then are they 
able to interpret the high Q feature as arising on account of 20 uniformly spaced 
sparks on the polar cap, rotating very stably round the magnetic axis, once every 
~ 37.35 pulsar periods (~ 41 s). Only after this are they are able to study the issue 
of aliasing and claim to have resolved it. For a brief and clear summary of the 
crucial a rgument of DR200 1 regar ding aliasing, see the second last para on page 
739 of Edwards & S tappers! (12002b . and the corresponding footnote. 

This is the first paper in the series that critically analyzes the DR1999 and 
DR2001 papers regarding the issue of aliasing. It argues that they are not justified 
in using this model as an essential input to resolve the issue of aliasing in PSR 
B0943+10; there are several factors within the framework of that model, that are 
capable of perturbing very steadily drifting sub pulses, thereby reducing the Q of 
the corresponding spectral feature to values much smaller than 500. The above 
factors are not important for the canonical drifting pulsars, in which the Q is not 
higher than « 10. In this paper the RSI 975 model also implies the subsequent 
modifica tions made to it by R uderman and colleagues, while retaining its essential 
features dCheng & Rudermanl[l977[ll980l: henceforth CR1977 and CR1980). 

RSI 975 derive the characteristic period P3 for a sub pulse position to recur as 



P 3 = — x 
N 



2nR s 1 . , 
x — periods, 

Vd P 
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where R s is the radial distance of the spark from the magnetic axis, N is the number 
of sparks laid out uniformly on the circumference 2nR s , P is the pulsar period, 
and Vd is the tangential speed of the spark, due to drift in crossed electric E and 
magnetic B fields in the vacuum polar gap (fig. [U; c is the speed of light. The 
frequency of the feature in the longitude resolved fluctuation spectrum is then 
1/P 3 [cpp]; this paper discusses the stability of this quantity. 

Let the parameters R s , Vj and N of eq. (1) vary independently with time, due to 
reasons described below, with standard deviations cr Rs , cr Vd , and cr N , respectively. 
Then the expected Q of the spectral feature is given by 

-1/2 



Q 



+ 



°~V d 

V d 



+ 



N 



(2) 



2 Radial drift of spark in RS1975 

First, strict RS1975 scenario is assumed; i.e., the vacuum gap breaks down due 
pair production in the gap B field, and not due to pair production in the Coulomb 
E field of Fe ions (CR1977). Then the velocity of radial shift (or spread) of the 
spark is given by (CR1980) 

v.. If 

where h is the height of the vacuum gap and S is the radius of curvature of the B 
field lines. 

For the required B field line curvature of S » 10 6 cm (RS1975, CR1977, 
CR1980), the time to cross the polar cap (« 10 4 cm) is « 100 micro seconds 
(jus); see CR1980 for details. Therefore in the duration of a spark, which RS1975 
estimate to be 1 [is to 10 /is, which also tallies with the micro structure time scale 
expected in a pulsar with period 1.0977 s (PSR B0943+10), one can expect the 
spark to have traveled « 1 meter (m) to 10 m. Thus one can expect cr Rs to be >> 1 
m due to this process. 

Now, DR1999 invoke 20 sparks, each of size « 20 m, uniformly spaced apart 
from each other by « 45 m, on a circumference of radius w 145 m on the polar 
cap. Thus the spark might have moved radially in wards by a fraction > 1 /20 of 
its size during its own development, leading to R s /o~r s < 145, implying that the Q 
of eq. (2) might be limited to the value < 145. 

CR1977 discuss this issue and propose two solutions for dealing with it, which 
is a departure from the strict RS 1 975 model. Their proposals probably do not work 
for high Q drifting as assumed in PSR B0943+10. 
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Figure 1: Similar to fig. 1 of CR1977. A high energy photon initiates a spark 
at position 1, by forming an e + - e~ pair due to the Erber process (see RSI 975), 
in a vacuum polar gap of height h. The e~ moves along the curved magnetic B 
field line, accelerated down wards by the electric field E, and emits a high energy 
photon due to curvature radiation. This travels tangentially to the original B field 
line and creates another e + -e~ pair at position 2. The e + of this pair creates another 
pair at position 3, and so on. Thus the spark shifts radially to wards the magnetic 
axis with speed V s , while simultaneously drifting tangentially with speed V d in the 
crossed E and B fields. R s is the radial distance of the spark from the magnetic 
axis (which is not shown here). 
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First they propose that the B field line geometry could be sufficiently compli- 
cated to stop such a radial shift; but this might also make the canonical tangential 
drift non-uniform, and once again decrease the Q of the spectral feature due to 
drifting. CR1977 also discuss another fundamental problem with this scenario. 

Second, they propose that the relativistic electrons that penetrate the surface 
of the neutron star might initiate a particle shower which might eventually lead to 
a high energy photon reentering the vacuum gap, to produce e + - e~ pairs. This 
allows the spark to be initiated on the original B field line in the vacuum gap, 
because one has now diluted the role of curvature radiation in the gap break down 
process. This might cause two problems for very stable (high Q) drifting: (a) 
unless this process is highly probable and occurs uniformly in space and time, it 
might cause erratic sparking in the vacuum gap, which might lead to a reduced Q 
of the drifting spectral feature, due to, say, random variations in shape, size and 
position of sub pulses; for example, small variations in the angle of entry into the 
gap of these high energy photons may change the position of the centroid of the 
spark; (b) unless the curvature related pair production ceases completely, one can 
not avoid altogether the radial shifting of the spark. This may lead to an effective 
widening of the spark area along the radial direction, if not to an outright shifting 
of the spark, which may also reduce the Q. 

CR1977 also propose an even more radical departure from the strict RS1975 
model. They propose that alternatively (to the Erber process) the gap can break 
down by e + —e~ pair production in the Coulomb field of Fe ions, which are released 
into the gap by the hot polar cap surface, by thermal x-ray photons from the same 
hot polar cap. They call this the Lorentz boosted pair production. The relativistic 
electrons returning to the polar cap surface maintain it hot. In this scenario also 
there are two problems. First, it is not clear what will be the final space distribution 
of the charges making up the spark - it may be much wider than « 20 m, and its 
centroid may differ from spark to spark by more than 1 m in position. Second, 
even if a narrow spark was drifting, and it carried underneath it the local hot spot, 
the cooling time for the polar cap surface is fairly long, of the order of w 0.2 s to 
even 10 s (CR1980). Therefore the spark may not drift as one observes, but may 
merely widen in the tangential direction. The non- uniformity of this widening 
may once again causing problems for high Q drifting spectral features. 



3 Change in drift speed of spark 
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So far the drift speed V d was assumed to be unvarying. However in the strict 
RSI 975 model, it may also depend upon the radial distance R s (see eq. (30) of 
RS 1975 and the subsequent discussion). Thus a perturbation in R s could also lead 
to a perturbation in Vd, further reducing the observable Q in eq. (2). 

Apart from the above, there appear to be very few perturbations that directly 
affect the drift speed Vd, at least to the levels discussed here. However, the curva- 
ture drift plasma process might set an upper limit to the observable Q. This arises 
due to the motion of charges along a curved pat h. The corresponding drift speed 



(in the tangential direction) is then given by (Landau & Lifshitz 1975; Chen 1984) 



V c = -T L U + Iv 2 l]\ SxB \- (4 > 



Scd b \ 11 2 

where y = 1/ a/1 - v 2 /c 2 , V\\ and v ± are the speeds of the charges parallel and 
perpendicular to the B field, co B = e\B\/m e c is the cyclotron frequency, and S and 
B are unit vectors normal and parallel to the B field line, respectively. Setting 
v,| * c, v ± « 0, oj b « 10 19 Hz, S « 10 6 cm, and y * 10 2 to 10 3 (RS1975, 
CR1977, CR1980), one obtains V c » 10~ 4 m/s to 10 3 m/s. Comparing this with 
the Vd = InRJAl « 22 m/s speed of sparks required by DR1999, the perturbation 
to Vd can be estimated to be o- Vd /V d of » 0.5 x 10~ 5 to 0.5 x 10~ 4 . Currently the 
latter number would imply a Q that is a comfortable factor of « 40 higher than 
the Q » 500 claimed by DR2001. So in principle it might appear that one can 
ignore this perturbation to Q in eq. (2). 

However, DR2001 claim that "longer fast Fourier transforms yield a progres- 
sively narrow feature" in their section 3; and they have used a data of length 
816 periods only. Even then they claim that the spectral feature "is only partly 
resolved". One should worry about the possibility that future (and longer) obser- 
vations might bring the observed Q of PSR B0943 + 10 uncomfortably close to the 
above upper limit. Of course, one should eventually do a full relativistic quantum 
calculation of eq. (4), and worry about the possibility that future theoretical inputs 
might also change the above upper limits. 

Furthermore, the above drift is charge dependent, i.e., the spark might get 
charge separated. This may not only widen the spark as it develops, it may also 
have very important consequences for the simultaneous existence of several iso- 
lated and independent sparks on the polar cap, as required by DR1999. 
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4 Change in number of sparks 



In eq. (2) a random change of number of sparks N by even +1 causes the Q to 
be « 20, much less than the >> 500 claimed by DR2001. Therefore any change 
in N is not allowed in the scenario of DR1999; which is why they claim that 
in PSR B0943+10 "20 sub beams rotate almost rigidly, maintaining their num- 
ber and spacing despite perturbations and that "the 20- fold pattern is always 
maintained strongly as a stable configuration" (on the right side of page 1010 in 
their paper). Now, there are two important elements of the physics of pulsars that, 
together, might not allow the above requirement of DR1999. 

First, the polar cap surface is supposed to be a good electrical co nductor; this is 
a nec essary requirement for the basic electrodynamics of the pulsar (iGoldreich & Julian 
1969). Even if the strong B field makes the transverse conductivity at the surface 
small, there must be a capacitive conductivity (i.e., conductivity with a time de- 
lay caused by, say, returning currents) near the polar cap surface, for the basic 
electrodynamics to hold true. 

Second, the gap break down process is supposed to be an extremely sensitive 
one. RS1975 state in their section 5 that "the extreme sensitivity of the discharge 
growth to small variations in the gap thickness may be reflected in the erratic 
nature of pulsar sub pulses". Later modifications to the RS1975 model dilute 
this sensitivity, since RS1975 assumed a gap potential difference that varies as h 2 , 
while space charge limited flow in the gap makes it vary only as h (eqs. (5) and 
(6) of CR1980). However, these later modifications also made the gap potential 
difference very sensitive to the thermionic emission from the polar cap surface, 
which depends exponentially upon the surface temperature, which in turn depends 
upon the return current heating the surface (CR1980). CR1980 state at the end of 
their section 3 that "This extreme sensitivity of the potential drop to the surface 
temperature makes the gap thermostatically regulating 

Given these two aspects of pulsar polar gap physics, the simultaneous exis- 
tence of « 20 stable sparks on the polar cap surface may not be assured. The 
spark growing the fastest on the polar cap might reduce the overall gap poten- 
tial difference, thereby inhibiting the growth of other sparks; at least one obser- 
vational paper has claimed to have seen competition in in tensity between drifting 
sub pulses (of PSR B0031-07: IVivekanand & Jo shi 1999). Observations certainly 
allow fo r 2 or 3 simultaneous spark s on the polar cap (see the earlier references; 
also see Vivekanand & Joshi (fl999k However, to the best of the knowledge of 
this author, no one has seen as many as >> 5 simultaneous sub pulses within a 
period. 
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In view of the above, the existence of « 20 simultaneous sparks on the polar 
cap is highly suspect, at least within the current RSI 975 framework. 

5 Thermostatic regulation of gaps 

CR1977 state in their abstract that "A potential drop above a polar cap may arise 
from (A) the inertia of freely released ions in a space-charge limited flow, or from 
(B) surface binding of ions leading to gap formation". CR1980 develop this idea, 
and claim that the gap potential will be, in general, intermediate between the value 
for a true vacuum gap and the value for a space charge limited flow. Such a 
situation can easily break into oscillations - a perturbation leading to an increased 
surface temperature will cause higher ionic current from the polar cap, resulting in 
a higher space charge and therefore a lower gap potential. This would reduce the 
(return) heating current onto the polar cap, thus cooling it, and reducing the ionic 
current. And vice versa when the perturbation decreases the surface temperature. 
From eq. (14) of CR1980, one can estimate the (cooling) time for the polar cap 
temperature to decrease by a factor of 10 to be ~ 0.2 s; from their eq. (16) 
the (heating) time for the polar cap temperature to increase by a factor of 10 is 
w 0.008 s. However CR1980 state that the actual times may be larger due to the 
approximations made. 

Therefore any spot on the polar cap, across which a spark has drifted, might 
have in general a quasi-periodic temperature variation, with a rise time of a few 
times « 0.008 s and a fall time of a few times w 0.2 s; the period of this variation 
will be the time for a spark to recur at the same spot, which is P 3 pulsar periods. 
How periodic or quasi-periodic will depend upon the actual period and the rise 
and fall time profiles, and their variation. The above time scales will change if 
the threshold is set not to increase or decrease of temperature by factor 10, but by 
some other factor. 

Now the gap potential is exponentially sensitive to this polar cap surface tem- 
perature variation (CR1980). So one may expect to see a corresponding quasi- 
periodic intensity modulation of the drifting sub pulse. In the spectrum of intensity 
fluctuations, this will be evident as side bands beating with the drifting spectral 
feature; their width and spacing will depend upon P 3 and Q of the modulating sig- 
nal. Indeed, DR1999 claim that in PSR B0943+10 the drifting feature at period 
P 3 S3 1.87 periods is amplitude modulated by a secondary periodicity of period 
« 37.35 pulsar periods (« 41.0 s). However, this amplitude modulating signal 
also appears to have a very high Q, of the order of 500 to 1000. Whatever be the 
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physical reason for such amplitude modulation, one should ensure theoretically 
that such an observation is possible; i.e., the quasi-periodic temperature varia- 
tion discussed above, which must be surely operative in the scenario of CR1980, 
must have a very high Q, at least higher than 500 to 1000. If not, then a low Q 
quasi-periodic temperature variation will broaden any narrow drifting feature in 
the fluctuation spectrum. 

From the theoretical viewpoint, a problem that the scenario of DR1999 might 
pose, is that any spot on the circumference 2nR s is heated 20 times faster than 
anticipated earlier, i.e., once every 1.87 pulsar periods. This might keep the av- 
erage polar cap temperature hotter than earlier anticipated; please note that the 
heating plus cooling timescale may be comparable to this periodicity. Then one 
might expect to see a significant amount of pulse nulling in the 816 periods of 
PSR B0943+10, whereas not even one pulse is observed to be nulled (DR2001). 



6 Conclusion 

This paper examined the issue within the frame work of the RS 1975 model, along 
with modifications made by CR1977 and CR1980. However, It is well known that 
the circular symmetry for the polar gap assumed in RSI 975 is an idealization. A 
realistic pulsar, having mis aligned rotation and the magnetic axis, has a north- 
south electric field on the polar cap dMichelll99ll) . This causes the sub pulse drift 
to be mainl y in the east-wes t direction on the polar cap, instead of being circular 
(fig. 2b oflRudermanl ll976h . This may modify significantly the conclusions of 



DR1999, where some sort of circular symmetry is essential. 

This paper is not a criticism of the RS 1975 model. It merely disagrees with its 
use, given its limitations, as an essential input to resolving the issue of aliasing of 
the high Q spectral feature in PSR B0943+10, even if that can actually be asso- 
ciated with the drifting sub pulse phenomenon. If such an association is actually 
invalid, then as to what actually is this spectral feature is a speculation beyond the 
scope of this work. 
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